Mon. Not. R. Astron. Soc. OOP. [TlfTTI C2008') Printed 2 Marcli 2009 (MN 1^1^ style file v2.2) 



Mapping the magnetosphere of PSR B 1055— 52 



o 
o 



6 



(N 
> 

oo 
(N 

(N 
O 

o 



Patrick Weltevrede^* and Geoff Wrigfit^'^ 

''-Australia Telescope National Facility, CSIRO, P.O. Box 76, Epping, NSW 1710, Australia. 
^Astronomy Centre, University of Sussex, Palmer, BNl 9QJ, UK 



ABSTRACT 

We present a geometric study of the radio and 7-ray pulsar B1055-52 based on recent 
observations at the Parkes radio telescope. We conclude that the pulsar's magnetic 
axis is inclined at an angle of 75° to its rotation axis and that both its radio main 
pulse and interpulse are emitted at the same height above their respective poles. This 
height is unlikely to be higher or much lower than 700 km, a typical value for radio 
pulsars. 

It is argued that the radio interpulse arises from emission formed on open fieldlines 
close to the magnetic axis which do not pass through the magnetosphere's null (zero- 
charge) surface. However the main pulse emission must originate from fieldlines lying 
well outside the polar cap boundary beyond the null surface, and farther away from 
the magnetic axis than those of the outergap region where the single 7-ray peak is 
generated. This casts doubt on the common assumption that all pulsars have closed, 
quiescent, corotating regions stretching to the light cylinder. 
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1 INTRODUCTION 

PSR B1055-52 is a middle-aged (0.5 Myr) energetic {E = 
3 X lO''* erg s^^) pulsar which has been detected as a 
source of pulsed 7-rays. In the ir study of energetic pulsars, 
IWehevrede fc JohnstonI l|2008l ) singled out PSR B1055-52 
for comment since it is exceptional in having both an inter- 
pulse and wide radio components, with the interpulse (IP) 
well separated from the main pulse (MP). In general the 
presence of an IP is a strong clue to a pulsar's geometry, and 
suggests that a study of the radio emission of PSR B1055- 
52 may give insight into the structure of the magnetosphere 
and the origin of its 7-rays. 

PSR B1055-52 was disco v ered w ith the Molongolo tele- 
scope by IVaughan fc Largel l|l972h . who noted a possi- 
ble IP in one of their records. The IP was confirmed by 
iMcCulloch et all (|l976l ') , who commented that the pulse 
profile resembles that of the Crab pulsar. iRankinI (Il983l ) 
noted that the MP-IP separation shows no variation with 
frequency, which is additional evidence that the pulsar's 
magn etic axis is h ighly inclined to the rotation axis. 
iLvne fc Manchcste3 l)l988l ) concluded that the observed po- 
larisation position angle (PA) swings of the MP and IP 
suggest an angle of inclination of the magnetic axis of 75°, 
without ruling out higher values around 90°. This general 
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view was supported in a complex study published bv iBigg j 
(1990a). 

PSR B 1055-52 is relatively close, which makes it pos- 
sible to observe its emission over a wide range of wave- 
lengths. The pulsar's DM distance is only 0.72 ± 0.2 kpc 
l|Cordes fc Lazioj 20021). roughlv c onsistent with the X-ray 
data (e.g. iQgelman fc Finlevlll993h and the non-thermal ra- 
dio source around PSR B1055-52 (jCombi et al.| [l997'). 

The EGRET detector (Energetic Gamma-Ray Exper- 
iment Telescope) on board of the Compton Gamma Ray 
Observatory (CGRQ) revealed that PSR B1055-52 is a 
pulsed 7-ray source (|Thompson et af] Il999l ) . Unlike most 
7-ray profiles (but in common with PSR B1706-44), that 
of PSR B 1055-52 does not show two narrow peaks sepa- 
rated by 0.4-0.5 in phase, but has a broad 7-ray profile with 
possible sub-peaks at the leading and trailing edge of the 
profi le which peak 0.25 and 0.05 in phase before the radio 
MP (|Thompson et al.lll999l ). The leading side of the profile 
has also been seen with the COMPTEL detector (0.75-30 
MeV), while only the trailing side of the profile is seen above 
2 GeV by EGRET, suggesting a longitude dependent spec- 
tral index. 

Another feature which makes PSR B 1055-52 an inter- 
esting source at high energies is its spin-down age of 0.5 Myr, 
making PSR B1055-52 (at the time of writing) the oldest of 
the known 7-ray pulsars. Its period of 0.197 seconds makes 
it also the slowest (excluding Geminga) . This status as being 
both an 'old' young pulsar and, at the same time a 'young' 
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Figure 1. The pulse profile of the 1369 MHz timing data of PSR B1055— 52 at 20 cm showing total intensity (solid line), linear polarization 
(dashed line) and circular polarization (dotted line). 



old pulsar (most radio pulsars are over 1 Myr old) suggests 
that it may possess characteristics of both and give insight 
into the link between high-energy production in the outer 
magnetosphere and radio emission above the polar cap. It 
has famously been dubbe d as one of the 'Three Musketeers' 
l|Becker fc Trumped [l997l ). along with PSR B0656+14 and 
Geminga, because these pulsars were found to emit solely 
thermal X-rays which are interpreted as the product of ini- 
tial neutron star cooling. 

The geometry linking t he high-energy and radio pro- 
files was recently studied bv IWang et al.l ( 2006 ) , who made 
use of the PA analysis of iLvne fc Manchester 1 19881 . They 
concluded that that the radio MP and IP are emitted from 
separate poles with the 7-ray emission coming from the same 
pole as the radio MP. They claim that the IP and MP come 
from different heights, but do not consider the effect of re- 
tardation and aberration on the PA-curve. In their model 
the radio emission occurs on fieldlines which are confined to 
a polar cap bounded by the last closed fieldlines which, in a 
dipole geometry, touch the light cylinder. Here we use new 
data to re-analyse the PA-swings and to create a geometry 
which is fully consistent with retardation and aberrational 
effects. A full understanding of the geometry is an important 
prelude to undertaking an analysis of the pulsar's complex 
single pulse behaviour, which will appear in a subsequent 
paper. 

The organisation of this paper is as follows: in Sect. 
[2] details about the radio data used are provided, in Sect. 
[3] the properties of the pulse profile of PSR B 1055-52 are 
described and its most likely geometry is derived in Sect.]!] 
In Sect. Owe identify the active fieldlines of the MP and IP 
and the implications of our results are discussed in Sect.|Sl 



2 OBSERVATIONS 

The pulse profile of PSR B1055-52 was obtained at three dif- 
ferent frequencies by the Parkes 64-m radio telescope located 
in Australia. During an observing session carried out from 



2006 August 24 to 27 data was recorded using the H-OH 
receiver (1369 MHz centre frequency, 256 MHz bandwidth 
split into 1024 freque ncy channels wit h an equivalent system 
flux density of 43 Jy; |johnstonll2003 ) and the 10/50 cm re- 
ceiver (3094/653 MHz centre frequency, 1024/64 MHz band- 
width split into 1024/512 frequency channels with an equiv- 
alent system flux density of 49/57 Jy on a cold sky). The 
signals from the two linear feeds of the receiver were con- 
verted into Stokes parameters, resampled and folded at the 
pulse period by a digital filterbank using 512 bins across 
the profile. In addition we also used the pulse profile with a 
slightly higher time resolution (1024 bins across the profile) 
obtained using the 20 cm multibeam receiver (1369 MHz 
centre frequency, 256 MHz bandwidth split into 1024 fre- 
quency channels with an equivalent system flux density of 
35 Jy on a cold sky) for the Fermi timing program. This pro- 
gram started in April 2007 with the aim to provide timing 
solutions which can be used to fold 7-r ay data obtained by 
the Fermi satellite (|Smith et al.ll2008l ). The proflle is the 
sum of 33 individual observations with a total length of 
6. 6 hours. For details about the data-processing we refer 
to I Weltevrede fc Johnstoiil l|2008l ). 



3 PULSE PROFILE 
3.1 Profile morphology 

The shape of the pulse proflle of PSR B1055-52 is quite 
comple x, as one can see in Fig. [1] In a similar way to iBiggsl 
(|l990al ) (who did not detect the first component of the MP), 
we call the four most prominent peaks of the MP peaks I 
to IV and the three peaks in the IP V to VII. Long integra- 
tions are required to obtain stable profiles, with evidence for 
persistent in tensity changes in component HI (also noted by 
lBiggs|[l990al ). The peak-to-peak separation between the MP 
and IP is ~ 159°, well below 180°, so it is not immediately 
obvious that we are dealing with emission from both poles 
of the neutron star. 
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Figure 2. Top panels: The pulse profiles of PSR B1055-52 at 
three observing frequencies drawn with different vertical offsets. 
FYom top to bottom the wavelength are f cm, 20 cm and 50 cm. 
Bottom panels: The PA of the linear polarization and an RVM 
fit. The offsets are the same as in the top panels. These are the 
August 2006 observations. 



As one can see in Fig. [T] the profile of the MP is highly 
linearly polarized. This is also the case for the leading edge 
of the IP, but at later puls e longitudes it become s almost 
completely depolarized (e.g. iMcCuUoch et al.l[l978l '). Notice 
that all three components of the IP have a distinct linear 
polarization peak, although they are much weaker for com- 
ponents VI and VII (Fig. [T]). The profiles of the degree of 
linear polarization at 10 cm and 50 cm (not shown) look 
very similar to that at 20 cm. 

Component II of the MP has a significant amount of 
circular polarization, although it is weak. There is a hint 
that the trailing part IP also has a significant amount of 
circular polarization, but with an opposite sign. 



3.2 Frequency evolution of the pulse profile 

The pulse profiles of PSR B1055-52 at different frequencies 
are shown in Fig. [51 The profiles are aligned and normalized 
using the peak of component II which is put at 90° pulse 
longitude. One can see that there is not much evidence for 
frequency evolution in the relative position of the compo- 
nents, but the peak amplitudes do evolve with frequency. 
In the MP the first component disappears at 50 cm while 
components III and IV become relatively stronger. The IP 
shows a similar frequency evolution with the leading half 
of the profile becoming weaker relative to the trailing com- 
ponent. In relative terms, the IP itself is strongest at high 
frequencies. While the components of the IP are sharper at 
10 cm this is less clear in the MP. Notice also that compo- 



nent VII is the only component which changes its centroid 
position with respect to the other components. 

There is also no evidence of the full widths of the MP 
and IP (or their separation) being frequency-dependent, po- 
tentially an important hint as to the pulsar's geometry. A 
profile a t 170 MHz, as cribed to McCuUoch et al. and pub- 
lished in iRankinlll 19831 ). shows a merging of components III 
and IV in the MP, but little overall widening of the profiles. 



4 GEOMETRIC IMPLICATIONS 

The PA-swing is plotted in the bottom panel of Fig.[2]for the 
three frequencies. We found a significant offset of the 50 cm 
data with respect to the other wavelengths wh en we use the 
rotat ion measure (KM) of 47.2 ± 0.8 rad m~^ (|Tavlor et all 
Il993l ). By applying an RM of 46 rad m~^ the PA-swings at 
the different frequencies line up, suggesting that the shape 
of the PA-swing, like that of the profile, is independent of 
frequency. This value of the RM is consistent with th e value 
of 44 ± 2 rad m"^ measured bv lNoutsos et al.l (|2008l ') . 

The shape of the PA-swing can be used to constrain 
the angles between the magnetic and rotation axis (a) and 
the angle between the line of sight and the rotation axis 
[C, = a + (3, where (3 is the impact parameter). According to 
the r otating vector model (RVM; iRadhakrishnan &: Cookg 
!l969!) the PA [ip) can be described as a function of the pulse 
longitude (0) by 

sin a sin (</) - (/jo) , , 

tan (?/) - i/>o) = ^— -r^. '—t- — (1) 

sm Q cos a — cos Q sm a cos (0 — (po) 

where ^po and (j^o are the PA and the pulse longitude corre- 
sponding to the inflection point of the PA-swing. Using this 
definition implies that we follow the observer's sign con- 
vention f or the PA , which differ s from that introduced by 
[Pa m o ur fc Tavlod fl992 ) and . Everett fc Weisberj (|200ll ). 

If the emission height r is zero, or at least negligible 
compared with the light cylinder radius Rlc, the inflec- 
tion point of the PA-swing coincides with the flducial point 
(defined as the pulse longitude where the magnetic pole at 
the surface is directed towards us). However, for non-zero 
emission heights this is no longer the case. As the emission 
height increases the total intensity appears to shift towards 
earlier pulse longitudes due to aberration (by r / Rlc) and 
the decreasing light travel distance to the observer (also by 
t/Rlc), while the PA-swing shifts to later pulse longitudes 
with respect to the fiducial point (by 2r/RLc)- Thus the 
total relative shift between the steepest gradient of the PA- 
swing with respect to the emission beam's centroid is given 
by 

A(I>BCW = 4,r/RLC (2) 

(|Blaskiewicz et al.lll99]] : lDvksll2008l ) , where r is assumed to 
be small compared with Rlc- 

In PSR B 1055-52 the degree of linear polarization is 
low in the trailing half of the IP, but fortunately it is 
enough to obtain the PA. This means we have a full PA 
coverage at virtually every pulse longitude where emission 
is detec ted at a far better t i me re solution than was avail- 
able to iLvne fc Manchester! l|l988l ) . This enables us to re- 
examine and place constraints on the geometric solution of 
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this pulsar. In this section we formally confirm that a nearly- 
orthogonal solution is preferable to a nearly-aligned solution 
and exploit the 7-ray profile to establish an integrated ge- 
ometry. 



4.1 Nearly-aligned rotator? 

We first consider the possibility that PSR Bf055-52 might 
be a nearly-aligned rotator. In such a geometry the MP and 
IP originate from the t wo edges of a hollow cone -like beam 
from a single pole (e.g. [Manchester fc Lvnj|l977h . The sep- 
aration between the two components can be large, but, as 
observed, will not necessarily be 180° in pulse longitude. 
The large widths of the components together with the steep 
trailing edge of the MP and the leading edge of IP are also 
typical for such geometry, making this model a promising 
possibility to consider. 

However, for a nearly-aligned rotator one expects the 
PA-swing to be flatter than is observed. In the RVM model 
the steepest gradient is associated with the line of sight pass- 
ing the flducial plane. In the case of an aligned rotator this 
plane should be at a pulse longitude between the MP and 
IP, hence the observed PA-swing will be relatively flat. As 
has been pointed out fbv IWang et al.l |2006| for example), 
the steepness of the PA-swing observed for both the MP 
and IP rules out the one pole geometry. It should be added 
that this conclusion is also valid even when 90° jumps in 
the PA-swing caused by different plasma modes dominating 
in diff erent parts of the pulse proflle (e.g. [Manchester et ahl 
1 19751 ) are considered. 

The PA-swing inflection point can in principle be shifted 
from its unshifted pulse longitude between the MP and IP 
(as would be expected for an aligned rotator) to appear at 
the pulse longitude of one of the pulse components (MP or 
IP) when the aberrational shift of Eq. [2]is considered. This 
effect could therefore in principle make the PA-swing steep 
in one of the two components. However, it cannot make them 
steep in both since the RVM model predicts a steep gradient 
at only one point in the PA-curve, thereby failing to explain 
the observed PA-swing. 



4.2 Nearly orthogonal rotator 

In searching for a self-consistent solution, we follow 



IWang et al.l l|2006l ) and consider a model where the proflle of 
PSR B 1055-52 is generated by emission from an orthogonal 

rotator. 

Although iMcCulloch et al.l (|l978f ) claimed that the 
PA-swing ca nnot be extrapola ted between the MP and 
IP, we (like Wang et al] I2OO6I with inferred data from 
iLvne fc Manchesterll 19881) had no problem in fltting the PA- 
curve with the RVM model (Eq. [ij. In fact, the reduced 
of 2.7 is quite good, although not perfect (Fig. The low- 
est is found for omp = 75° and (3mp ~ 36.1° (where 
aMP and Pmp are the MP values of a and /3), suggesting 
that the orientation of the magne tic axis is ind eed close to 
orthogonal and in agreement with IWang et al.1 (2006). This 
implies aip = 180° — qmp = 105° and an impact parameter 
at the IP of /3iP = /3mp + 2aMP - 180° = 6.1°. Note that 
these fits require that all the emission at each pole comes 
from the same height, a point to which we will return in 




90 

a (deg) 

Figure 3. The reduced grid obtained by fitting the RVM to 
the PA-swing of both the MP and IP of the 1369 MHz timing 
data. Black corresponds to 2.7 and white to 20. The cross indi- 
cates the position of the minimum of the surface (omp = 75° 
and /3mp = 36.1°). Solutions to the right of the right hand side 
diagonal solid line are excluded because they require a /3jp > 90° 
and solutions to the left of the left hand side diagonal solid line 
are excluded because they require a negative /3ip (corresponding 
to a declining PA-swing for the IP). Solutions below the dashed 
line (for q^mp < 90°) and above the dashed line (for omp > 90°) 
are excluded as they do not allow the 7-ray emission to be emitted 
from above the null line. 



Sect. 15.21 The best flt is shown in Fig. [2] as the solid curve. 
The biggest discrepancy between the flt and the data occurs 
in the trailing component of the IP, where the flt is too flat 
compared with the model. The steepest gradient in the MP 
{(j)o of Eq. [l| is found exactly at the location of the highest 
peak in the pulse proflle {(j> = 90° in the plots). 

The fact that the PA-swing can be extrapolated be- 
tween the MP and IP further implies that the emis- 
sion heights at both poles cannot be too large. If there 
were a large differential emission height, the PA-swings 
of the MP and IP would be shifted by different amounts 
jBlaskiewicz et al.l Il99ll : iHibschman fc AronsI I2OOII : iDvksl 
l2008l ). making it impossible to flt the PA-swing with a sin- 
gle RVM curve. By modelling these relative shifts we find 
that the lowest is found when the emission height of 
the IP is higher than that of the MP by 0.07iZLC (~ 700 
km) with an uncertainty of 0.10i?LC (~ 1000 km), which is 
consistent with the emission height being identical at both 
poles. The uncertainty is derived by finding the differential 
emission height at which the reduced of Ihe PA fit ha s 
doubled. Although the solution found bv lWang et al.1 (|2006D 
(who find an IP/MP emission height differential of between 
O.lORi^c and 0.20i?Lc) is in the allowed range, we will argue 
later that this solution is very unlikely. We therefore seek a 
solution with identical emission heights for the MP and IP. 



4.3 Implications of the 7-ray profile 

As lLvne fc Manchester! l| 19881 ) pointed out, and as is evident 
in Fig. |3l there are RVM solutions over a wide range of 
qmp values which have a comparable x^- However, we can 
utilize additional constraints to confirm that we have the 
right solution. By assuming that f3ip must be smaller than 
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90° (or else the full opening angle of the radio beam would be 
larger than 180°) and that /3ip must be positive (or else the 
PA-swing of the IP would decrease) only solutions between 
the two diagonal solid lines in Fig. [3] are allowed. 

A further and critical constraint can c ome from the ob- 
serve d phase and width of the 7-ray profile (|Thompson et all 
1 19991 ). Different magnetospheric locations are proposed for 
the 7-ray emission of pulsars. Early models put the acceler- 
ation zones at l ow altitudes above the magn etic poles (polar 
cap models; e.g. lDaughertv fc HardirL3ll996l ). More recently, 
it has been argued that the emission height could be consid- 
erably larger for fieldlines close to the l ast open fieldlines giv- 
ing ris e to the "slot gap" model (e.g. iMuslimov fc Hardind 
2004al ) and the "two-pole caustic" model ( Dvks fc RudakI 



3). A separa te class of rnodels are the so-calle d "out- 

ergap" models (|Cheng et all 1 19861 : iRomanil 1 19961), which 
place the acceleration zones above the "null lines" l|Hollowavl 
1 19731 ). These are the lines on which th e charge density of 
a pu lsar magnetosphere changes sign (iGoldreich fc JulianI 
1 19691 ). 

Numerical simulations show that only an outer gap 
model can reproduce the observed shape of the 7-profile 
of PS R B1055-52 to a good approximation ( Watters et al.l 
I2OO8I ). The two-pole caustic model, which has 7-ray pro- 
duction below the null line, predicts in general profiles with 
observed emission from both magnetic poles, which makes it 
difficult to reproduce the observed 7-ray pulse width. In or- 
der to reproduce the observed single-peaked 7-ray profile the 
two-pole caustic model predicts either a or C, values smaller 
than allowed by the PA-swing parameters between the two 
solid diagonal lines of Fig. |3l Polar cap models, which have 
low emission heights, also require too small a and values in 
order to get a broad enough 7-ray profile. On the other hand 
outer gap models predict, by definition, that only emission 
from one pole can be seen for any given observer, thereby 
qualitatively reproducing the observed 7-ray profile for a 
wide range of geometries. 

In order to reproduce the observed width of the 7-ray 
profile the gap thickness w (the fraction of the angle from 
the last closed fieldline to the magnetic axis) should be of 
the order of 0.18 (assuming that the 7-ray efficiency is equal 
to ui), showing that the gap is very large in comparison with 
most 7-ray pulsars l| Watters et al.l 120081 ). Large ou ter gaps 
are ty pical for highly efficient old 7-ray pulsars (e.g. iRomanil 
Il996l ). The outer gap model further predicts that two caus- 
tics are formed at the edges of the 7-ray beam, which could 
account for the sub-structure visible in the EGRET data. 

If we accept that the 7-rays are produced above the null 
surface of the MP (i.e. they are produced by the outer gap) 
large parts of the omp — /3mp plane can be ruled out. To- 
gether with the 0° < /3/p < 90° restriction we can eliminate 
the entire hashed regions in Fig. |3] from consideration. As 
one can see in Fig. [J] using these constrains the geometry of 
PSR B1055-52 is already constrained even without fitting 
the shape of the PA-swing of the IP as there is only a rela- 
tively small wedge shaped region of allowed solutions in the 
plot. This is important, because it shows that our conclu- 
sions do not depend strongly on the shape of the PA-swing 
of the IP which deviates slightly from the RVM prediction. 
We will therefore adopt the RVM solution with the lowest 
in this paper which has omp = 75° and /3mp ~ 36.1°. 
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Figure 4. The reduced grid obtained by fitting only the PA of 
the 1369 MHz timing data of the MP region. Black corresponds 
to 1.1 and white to 3.0. For an explanation of the lines see Fig.|3] 



4.4 The geometric picture 

The emission geometry in the fiducial plane is illustrated in 
Fig. [5l One can see that the line of sight passes the mag- 
netic axis of the IP much more closely than to the pole of 
the MP, consistent with the steeper slope of the PA-swing 
of the IP. Thus, the IP emission must arise from open field- 
lines as is generally thought to be the case for radio pulsars. 
However, the incidence of the line of sight with the magnetic 
axis of the MP is very large, which immediately leads to one 
of two possibilities. If we require the emission to arise from 
open fieldlines, then it must originate at a considerable frac- 
tion of the light cylinder. Otherwise we are forced to accept 
emission from presumed closed fieldlines at low altitudes. 

Fig. [S] also shows the null lines (dashed). These are the 
f2 ■ B = lines on which the charge density of the mag- 
netosphere changes sign. The tangent of the fieldlines are 
perpendicular to the rotation axis at the point where they 
intersect the null line. It immediately follows that for the 
assumed geometry we 'see' fieldlines above the null line for 
the MP (Fig. [S| , which suggests that the radio MP is pro- 
duced in a region where an outergap is likely to exist. It can 
be shown that this is a general rule. In any pulsar where 
emission from both poles are detected, at one of the poles 
our line of sight must see emission from the region above 
the null line, i.e. the likely site of any "outergap". Assuming 
that this outergap is the site of the observed 7-ray emission, 
it is no surprise that the phase range of the 7-ray profile 
overlaps that of the MP rather than the IP. We note that 
were we to choose an RVM solution with qmp > 90°, then 
the IP would be the pole which views the outergap region, 
and the 7-ray profile would be at the wrong phase. 

Most other known 7-ray pulsars are, unlike PSR B1055- 
52, typically found to have a single peaked radio profile with 
a large offset from the 7-ray profile. The lines of sight of 
these pulsars are also believed to have a large /3 value for one 
of the poles. This causes the line of sight to miss the radio 
beam of the magnetic pole for which the 7-rays are observed, 
hence a large offset is seen between the radio and the 7-rays. 
PSR B1055-52 is different because, for some reason, it has 
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Figure 5. The geometry of the dipole field of PSR B1055-52 
assuming omp = 75°. The rotation axis points upwards and the 
fight cylinder is indicated by the vertical lines. The line of sight 
(dotted line "LOS") makes a large angle (36.1°) with the mag- 
netic axis of the MP (right hand side of the plot), while this angle 
is much smaller (6.1°) for the IP (left hand side of the plot). The 
dashed lines indicate the location of the null surface, and the last 
closed fieldlines which conventionally define a corotating region 
are shown. At the MP we must be viewing fieldlines which lie 
above the null surface (in the shaded region) - and which will 
be well within the corotating region if the radio emission is, as 
usually assumed, close to the surface. The 7-ray profile, which is 
observed to overlap the MP radio emission, is most likely formed 
in the shaded area, which indicates the volume above the null 
surface where an outer gap could potentially be formed. 



a large enough radio beam to be seen by a line of sight with 
a large angle /3. 

One further conclusion can also be drawn. With /3mp 
being so large at the MP (36.1°), the corresponding Pip at 
the IP will be positive (since /3ip — Pmp + 2aMP — 180° = 
6.1°). This implies that at the IP we view the pole on the side 
away from the null line, and therefore do not see fieldlines 
which thread the outergap region. 



5 BEAM ILLUMINATION 
5.1 Beam mapping 

We now use the widths of the MP and IP pulse profile to 
place additional constraints on the magnetosphere's geome- 
try. The first point that needs to be considered is the devia- 
tion of the MP-IP separation from 180° in pulse longitude. 
This could in principle be caused by a differential emission 
height, but, as we argue in Sect. 15.21 this solution is very 
unlikely. This leaves us with the less satisfactory conclusion 



that at least one of the radio beams is not fully active (e .g. 
iLyne fc ManchestejflQSj : iKarastergiou fc JohnstonllioOTi '). 

Our next aim will be to investigate which parts of the 
MP and IP beam should be active in order to explain the 
observations. This is done by mapping the observed pulse 
longitudes onto the fieldlines which are pointing in our di- 
rection. The pulse longitudes of the edges of the MP and 
IP are therefore an important factor in determining which 
fieldlines are active. We take the positions of the steepest 
gradient of the RVM fit ((^0 and 180° + 0o, see Eq. as 
the reference. These are located in component II of the MP 
and at the trailing edge of the IP. The emission of the MP is 
located at a pulse longitude range of A0mp ~ —15° to +22° 
with respect to the first reference point and the emission of 
the IP confined to A(/)ip = —42° to +10° using the second. 

This can be used to identify the fieldlines which are 
active after correcting for the relative shift between the PA- 
swing (which is used as reference) and the total intensity 
for a non-zero emission height (Eq. [2} . This means the true 
longitude of the emission point with respect to the fiducial 
point in the corotating frame is A0 + A0bcw. If p is the 
opening angle of the fieldline at the emission height and Lp is 
the azimuthal angle of the fieldline with respect to the mag- 
netic axis such that zero corresponds to the fiducial plane, 
these are given by 



cos p — cos a cos ^ + sin a sin cos ( A^ + A(j!)bcw ) , 
sin ( sin {Acp + A^bcw ) 



sin (p = 



sm p 



(3) 
(4) 



where A0 — (p — (po is the observed pulse longitude differ- 
ence with respect to the reference point as dis cussed above 
l|Gil et al.lll983 : iGupta et al.l 12004 lDvksll2008l see also Ap- 
pendix |A] for a derivation). 

Note that p is only identical for the leading and trail- 
ing edge of a beam when the beam is symmetric about 
the fiducial plane. In that case A<;6ieading + A^bcw ~ 
— (A^traiiing + A^BCw) ~ \W , where W is the full pulse 
width. 

Assuming a dipolar field, the opening angle of the field- 
line p can be expressed in term s of the polar angle 9 (colati - 
tude) of the point of emission (|Gangadhara fc Guptall200ll ) 



tant 



2 tan p 



+ 1 2 + 



2 tan p 



(5) 



which then can be expressed in terms of the footprint pa- 
rameter s on the polar cap 



r 



(6) 



The magnetic axis corresponds to s = 0, and s = 1 corre- 
sponds to last open fieldlines. As will be discussed in Sect. 
5.31 the polar cap is not necessarily strictly circular if the 
magnetic axis makes an angle with the rotation axis. 

The geometrical solution for the leading and trailing 
edges of the MP and IP are shown in Fig. |6] as a function 
of the emission height (because Ai^bcw is a function of the 
emission height). These plots identify which part of the po- 
lar cap should be active in order to explain the component 
widths and separations in combination with the shape and 
shift of the PA-swing. One can see that a becomes very large 
when the emission height is very low. This is because there 
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Figure 6. The geometrical solution for which ficldlincs should be active in order to fit the PA-swing and the pulse profiles as function of 
the emission height. The fieldline is characterized by the footprint parameter s and the azimuthal angle ip. The solution for the MP are 
the thick curves and that of the IP are the thinner curves (assuming Qmp = 75°). The solid lines indicate the solution for the leading 
edge of the component and the dashed line that of the trailing edge. The dotted line indicates r = 0.07iJLC- 



is a trade-ofF between the emission height and the s value 
in order to get the required observed component width. At 
large emission heights there is a turnover and s increases 
with emission height. This can be understood in terms of 
the shift of the PA-swing with respect to the emission pro- 
file. In order to prevent the PA-swing from having its steep- 
est gradient far outside the beam, s has to increase. At the 
same time the fraction of the beam which is active has to de- 
crease in order to avoid generating profile components which 
are larger than observed. This can be seen in the ip curves, 
which approach each other for large emission heights such 
that only the trailing edges of the beams are active. 

5.2 The emission height 

It is not possible to determine the emission height based 
just on the shape of the pulse profile in combination with 
the shape of the PA-swing. Additional assump tions are re- 
quired to constrain the emission height further. IWang et al.l 
l|2006t ) make two assumptions in order to derive their solu- 
tion. Firstly, both the MP and IP beam are assumed to be 
filled (i.e. that Lp of the leading and trailing edge have the 
same magnitude but opposite sign) . This choice is not avail- 
able to us because we have insisted on a solution fully con- 
sistent with retardat ion and aberrational effects. The second 
assumption made bv lWanE~et al. (200(|) is that the s param- 
eteiQ of the MP is assumed to be 1. One can see in Fig. [6] 
that the trailing edge (dashed line) of the IP and especially 
the MP have to be emitted from fieldlines further away than 
s = 1 from the magnetic axis (which is true for all possible 
emission heights). This show s again that the solution pre- 
sented by I Wang et al.l l|2006l ) is not self-consistent. 

Fig.[6]shows that the solid and dashed s curves of the IP 
cross at r = 0.07J?lc- This corresponds to a beam which is 
completely filled such that the distances of the footprints of 
the active fieldlines to the magnetic axis are the same for the 
leading and trailing edge. This is consistent with the (p values 
at this height which have the same magnitude but an oppo- 
site sign. This happens at an s value slightly larger than 1, 

^ The s parameter corresponds to rj in I Wang et al. 1 l l200d) . 



indicating that the beam is slightly bigger than the conven- 
tional polar cap size. The same emission height corresponds 
to a special geometry for the MP as the solid ip line crosses 
zero at that height. This means that the MP beam is exactly 
half filled for the emission height at which the IP is exactly 
completely filled. This arises because at this emission height 

^loading + A(j>BCW = - (^trailing + A^BCw) for the IP. Fot 

the same emission height at the MP ^loading + A</)bcw ~ 0. 
There is no obvious reason why an exactly half filled MP 
should have physical significance, so it is presumably a co- 
incidence this solution occurs at the same emission height 
as the fully filled IP solution. Moreover, we cannot know if 
the emission height is exactly r = 0.07-Rlc- 

The trailing edge of the radio MP and the leading edge 
of the radio IP are relatively steep (Fig. [TJ . This suggests 
that we only see the trailing s ide of the radio MP and o nly 
the leading side of the IP (e.g. lLvne fc Manchesteilll988l ). If 
we accept this idea, then the emission height must be smaller 
than 0.07/?LC = 700 km, because for higher emission heights 
the IP's trailing edge (rather than leading) would corre- 
spond to an edge of the beam (see Fig.[S)). This would imply 
that the radio emission height of PSR B 1055- 52 is not ab- 
normal compared with other pulsars (e.g. .Blaskiewicz et al.l 
19911: iMitra fc RankinI l2002l : IWeltevrede fc JohnstonI l2008l : 
Kiiak fc Gillll997l ). 



If there is a differential emission height between the two 
poles, then that of the IP is larger (see Sect. 14.2")) . However, 
the IP's emission height cannot be larger than 0.07iiLC, or 
else its leading edge would no longer be the edge of the beam. 
The emission height of the MP could be slightly lower, al- 
though not too much as it would make its s value unrealis- 
tically large. Therefore roughly similar emission heights at 
0.077?LC (or slightly lower) is the most likely solution. 



5.3 Polar map 

Having chosen a preferred (maximum) height, we may use 
the geometric solution of Fig. [6] to generate a map of the 
polar cap of PSR B1055-52. The result is shown in Fig. [7] 
and combines the IP and MP on a single pole, thereby as- 
suming that the emission structure at both poles is identical. 
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Figure 7. The polar cap map of PSR B1055— 52 based on the 
assumption that both poles have identical structure so that we 
can "reflect" the observed IP mapping in the horizontal axis. The 
conventional polar cap size of an aligned rotator is indicated by 
the circle with radius s = 1. The dashed ellipse shows the polar 
cap including meridional compression. The fieldlines inside the 
kidney-shaped area are open and do not pass through any null 
line. The top and bottom solid curves are the lines of sight for the 
MP and IP respectively, where the active portions (assuming an 
emission height of 0.07-Rlc) the thicker parts of the curves. 
The location of the peaks are indicated as well (see Fig. [!}. The 
dotted semi-circle indicates approximately the footprints of the 
fieldlines on which 7-rays are thought to be produced. 



Obviously, because we only have information about the line 
of sight which passes over different parts of the polar cap 
of the MP and IP, we cannot know whether this assump- 
tion is true. It must be borne in mind that both the IP 
and MP are observed on the same meridional planes of ip 
(i.e. they both are viewed below the magnetic axis in Fig. 
[5|, and in Fig. [7] the IP line of sight traverse has been "re- 
flected" in the horizontal axis. This reflection ensures that 
a point in the polar map is associated with fleldlines which 
have identical geometric properties for both the MP and 
IP. The circle s = 1 corresponds to the footprints of the 
last closed fieldlines of an aligned dipole touching the light 
cylinder. For non-aligned dipoles the po lar cap shap e gets 
distorted by meridional compression (e.g. lBig3l990bl ). The 
predicted shape for a dipole field inclined at 75° is shown as 
the dashed oval in Fig. [T] However, it must be noted that 
even models which take into account the distortion of the 
magnetic field lines close to the light cylinder predi ct roughly 
circular beams (e.g. iRomani fc Yadigaroglulll995h . 

In Fig. [7] we also indicate the kidney-shaped region 
within which fleldlines of a dipole inclined at 75° do not 
cross the null surface inside the light cylinder. This shape 
is numerically derived by determining for each dipole field- 
line whether the null surface f2 ■ B = is penetrated within 
the light cylinde r. It can also be determined analytically 
l|Qiao et al.ll200j ). In the upper part of the polar map the 



null surface will lie at a great distance from the polar region 
(see Fig. [5]), but in the lower section it may be close enough 
to the star to set up an active outergap. We also indicate the 
approxim ation of the footprint s of the = 0.18 fieldlines, 
argued bv lWatters et al] l|2008l ') to be the fieldlines on which 
the 7-ray emission is produced (dotted semi-circle in Fig.[7|. 

Note that the locus of active fieldlines on which the IP is 
generated is symmetric about = and barely extends out- 
side s = 1. Furthermore, the IP fleldlines, apart from their 
outer edges, are contained within the kidney-shaped locus of 
the open flel dlines which do not cross a null surface. It was 
suggested bv lGoldreich fc JulianI (1969) that such fieldlines 
carry the outflowing current from the pulsar. 

By contrast, the footprints of the active MP fieldlines 
lie well outside s — 1, extending to radii more that double 
this value. Strangely, they cover only the latter half of a 
symmetrically-defined full beam. However, as pointed out 
above, the MP's fieldlines will thread the null surface close 
to the star's surface (at about 160 km) and the emission will 
occur in a plasma region usually considered to be inactive 
and corotating, so we can be sure that the MP is not a 
traditional polar cap profile. 

The solution shown is for r = 0.07J?lc, which is argued 
to be the maximum allowed emission height (Sect. [F!2)) . For 
smaller emission heights both the line of sight of the MP 
and the IP will move to larger values of s (see Fig. |6} and 
the active parts of the lines of sight will start leftwards of 
their current location (more negative values of ip). 



6 DISCUSSION 

In our study of PSR B 1055-52 we have focused on inves- 
tigating the geometry of its magnetosphere. Its wide MP 
and IP, together with a clearly-positioned single gamma- 
ray peak, place considerable constraints on any potential 
model. In an earlie r attempt to locate the emission regions, 
IWang et al] i200di) fitted an RV M model to published fig- 
ures i Lvne fc Manchester! Il988l ) for the pulsar's PA-swing 
and did not consider the constraints that aberration and re- 
tardation might place on the choice of emission heights. We 
have been able to exploit new data, gathered at Parkes in 
recent years, to verify the RVM fit and to attempt a self- 
consistent model. 

One of our principal conclusions is that both the MP 
and IP are emitted at roughly the same height. The optimum 
fit (qmp = 75°, Pmp = 36.1°) is in line with earlier results. 

This PA fit has a number of interesting implications 
for PSR B 1055-52 which are independent of assumptions 
about the emission heights. Since ^mp (i.e. aup + Pmp), at 
111° is considerably greater than 90°, our line of sight at 
the MP must see emission on fleldlines at locations which 
are above the "null" line (Fig. [5]), in or beyond territory 
normally ascribed to an outergap (as was also noted by 
iHirotani fc Shibatal |2003 '). In fact, on elementary geomet- 
ric grounds, it must be true for every pulsar with a detected 
IP that at one pole we are viewing above the null line and 
at the other below it — a point not widely realised before. 

However, in the case of PSR B1055-52 (but not neces- 
sarily for all pulsars with IP's), ^ip > aip and so our line 
of sight at the IP must cross the magnetic axis on fieldlines 
which do not thread the outergap at the opposite pole (Fig. 
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[7|. On these grounds alone, we can expect that emission at 
the MP and IP will reflect the different conditions of these 
two regions. This can be contrasted with the emission from 
PSR B1702-19, where at the MP, in one of its two possible 
solut ions, q:mp ~ 99° and C,mp ~ 91.5° (|Weltevrede et al.] 
I2OO7I ). Then, at one pole (IP) we must be viewing fieldlines 
above the null line and at the other (MP) equivalent field- 
lines before they cross the null line, with a greater likelihood 
of coordinated emission. 

Although the radio emission at both poles must come 
from about the same height, the precise height is not easy to 
determine. To first-order, aberration at a fixed height sim- 
ply shifts the entire PA curve to later pulse longitudes by 
an amount proportional to the height. Thus, as long as the 
emission height is small compared to the light cylinder ra- 
dius and unchanged across the profile, we have no way of 
knowing from the shape of the PA curve by how much it 
has been shifted. The exceptionally wide MP and IP lead 
us to believe that, whatever the assumed height, it is im- 
possible to confine the footprints of the emitting regions to 
a polar cap of conventional dimensions (i.e. defined by last 
closed dipole fieldlines touching the light cylinder). This re- 
quirement is particularly acute for the MP, whose trailing 
components must lie on fieldlines located about twice the 
conventional polar radius from the pole. Furthermore, since 
the fiducial point of the MP profile must lie ahead of the 
steepest gradient of the PA-swing, the observed MP width 
will never occupy much more than half its notional beam. 
Choosing a height less than about 500 km would force the 
polar cap size to escalate rapidly (see Fig. [Sjl and therefore 
a balance has to be struck. Here a height 7 per-cent of the 
light cylinder radius is assumed, so that the IP emission fills 
its beam and has footprints at around s = 1 on the oppo- 
site side of the pole (assuming identical polar caps at both 
poles). This choice also has the advantage that the steep 
edges of the MP (trailing) and IP (leading) correspond to 
beam edges (see Figs. |6] and [T]). We therefore conclude that 
the emission height is unlikely to be higher or much lower 
than 700 km, a typical value for pulsars. 

The kidney shaped central area of Fig. [T] which almost 
completely encloses the footprints of the IP's line of sight, 
highlights a significant consequence of the inferred geometry: 
the active fieldlines of the IP are open fieldlines which not 
only do not thread the outergap region beyond the null line 
close to the surface, but they also do not thread the null 
surface at any point (assuming a pure dipole geometry as 
in Fig. [Sjl before they leave the light cylinder. This suggests 
that along these fieldlines the pulsar may expel or accrete 
charged particles without having to encounter any charge 
sign boundary. Such fieldlines appear to bunch around the 
magnetic pole bounded by the radial parameter s — 1, the 
conventional limit for a polar cap defined by the last closed 
fieldlines touching the light cylinder. Thus the IP, if observed 
on a line of sight which missed the MP, would not suggest 
an unusual geometry. 

By contrast, the fieldlines generating the MP have polar 
footprints well beyond s = 1 (Fig. [7| , although the observed 
emission region is itself well above the null line (Fig. [5]). We 
are therefore led to the conclusion that in PSR B1055-52's 
MP we are seeing emission from regions of the magneto- 
sphere normally regarded as a closed and corotating. A polar 
cap radius of s = 2 is the footprint of a fieldline which closes 




rotation 
axis 



Figure 8. A schematic picture of a physical model for 
PSR B1055-52, deduced from its radio and 7-ray profiles. The 
corotation region, bounded by the fieldline s = 2, is shaded and 
lies well within region defined by the last closed fieldline to touch 
the light cylinder. It is surrounded by an "auroral" zone, which 
may or may not be entirely closed, but which must support elec- 
tric fields parallel to its magnetic fieldlines. It makes two dif- 
ferent intersections with the null surface, one very close to the 
star and a second no longer as distant from the surface as in a 
conventionally-defined corotation region. This may give rise to 
a complex outergap structure, with pair-created particles con- 
tinually reflected between the null surfaces. The 7-ray emission, 
appropriate to a "middle-aged" pulsar, is located above the out- 
ergap, stretching from the null line to the light cylinder. 



at a distance of only a quarter of the light cylinder radius 
and which intersects the null line at a height of only 160 km 
above the polar cap (Fig.jS)). In Fig.[7]we see that radio emis- 
sion occurs on fieldlines with parameters s — 1.5 (or smaller) 
to s = 2, indicating emission on fieldlines for a significant re- 
gion between the s = 2 surface shaded in Fig. |5] and the last 
closed fieldlines (s ~ 1). Radio emission from these fieldlines 
implies particle acceleration and hence the presence of elec- 
tric fields parallel to the fieldlines. So it is conceivable that 
such fieldlines are open and cross the light cylinder (as they 
do, for example, in the models of aligned magnetospheres 
with reduced corotation regions given by iGoodwin et al.l 
|2004 ). The closeness of the outergap to the neutron star 
surface may lead to an active outergap accelerato r merging 
with the polar gap (see iHirotani fc Shibatall2002l: iHirotanil 
l2006l : lMuslimov Hardinell2004bl : lQiao et al.ll2004l ') and the 
much reduced corotating region (shaded region in Fig. [SJ . 

Alternatively, the quasi-closed region between fieldlines 
s = 1 and s = 2 can be seen as an auroral zone, analogous to 
those known to be present in the planets, in which accretion 
from a d ebris disk m ay occur and modulate the radio signal 
(Cordes fc Shannon , 2008: L uo fc Melrose 2007 ). In fact, ac- 
cretion could modulate the entire outergap regio n through 
an injection of charge (jHirotani fc Shibatall20o3 l. Since all 
the surface s = 2 is now so close to the neutron star, it is 
not impossible to envisage interaction between the two rmll 
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surfaces which intersect it in the meridional plane (Fig. |8]). 
Possibly "gaps" arise at both at which pair-creation may be 
stimulated by thermal X-rays from the star's surface lead- 
ing to secondary particles "bouncing" between the gaps and 
stimulating the observed MP emission. 

Whatever conditions are present in the envelope be- 
tween s = 2 and s = 1, PSR B1055-52's 7-rays are al- 
most certainly generated relatively close to the light cylin- 
der above an outer gap which stretches from the last closed 
fieldline (whether s = 1 or s = 2) towards the light cylin- 
der. Unlike most other 7-ray pulsars, B1055-52 ha s a single 
broad peak which is shown by Watte rs et al.] (|2008l ) to match 
cl osely to the geo metric predictions of the outergap model 
ofE Oman ] (|l996h . lying on a fieldline with s = 0.82. Out- 
ergaps with 7-ray emission regions lying significantly above 
the last closed fieldline are thought to be a feature of rela- 
tively "old" pulsars, as PSR B1055-52 indeed is by the stan- 
dards of known 7-ray pulsars. The fact that the 7-ray pulse 
overlaps the phase of the MP indicates that both regions 
point to us at about the same time, which lends support to 
our model (Fig. [Sjl . 

Our conclusions about the geometry of PSR B1055- 
52 result from a rigid application of the RVM model to 
high-quality data, taking into account the expected ef- 
fects from aberration. We are convinced that this pul- 
sar emits radiation from fieldlines lying well outside the 
conventionally-defined polar cap boundary. Consequently, 
we would question the common assumption that all pul- 
sars have closed, quiescent, corotating regions stretching 
to the light cylinder. Possibly this only occurs after the 
outergap has become inactive, so that energetic (high E) 
pulsars can often have wide radio profiles, as observed 
I Weltevrede fc Johnstonll2008l ). Large polar caps might also 
help explaining the inconsistency found in calculating emis- 
sion heig hts of energetic pulsars using last-closed-fieldline 
methods l| Weltevrede fc Johnstonll2008l ). 

An unsatisfactory aspect of our model is the one- 
sidedness of the MP emission region. This might arise from 
the very different azimuthal trajectories which particles on 
the trailing and leading halves must take. On the trailing 
half they flow counter to the sense of the pulsar's rotation 
and can smoothly cross the light cylinder to the wave zone. 
However, those leaving the leading half of the polar cap will 
increase speed and hence Lorentz factor as they gain height 
on corotating fieldlines. It has been suggested that this may 
lead to additional pair-creation and high-energy emission as 
the parti cles free themselves from fieldlines near the light 
cylinder (jMestel et al.lll976l : Ida Costa fc KahnI Il982l ). but 
subsequent screening and its consequences for radio emis- 
sion have not been considered. 
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APPENDIX A: TRIGONOMETRIC RELATIONS 

We are interested to find the relation between the observed 
pulse longitude and the magnetic coordinates p (the angular 
distance from the magnetic pole) and Lp (the azimuthal an- 
gle around the magnetic pole). This coordinate transforma- 
tion can be derived using spherical geometry on the celestial 
sphere of the pulsar. Two trigonometric relations wh ich can 
be us ed are the sine and cosine rules for sides (e.g. ISmartI 
Il96d ) 

sin a sin b sin c , , , , 

(Al) 



sin A 



sin B sin C 
= cos b cos c + sin 6 sin c cos A. 



(A2) 

See the top panel of Fig. lAll for the definitions of the angles 
in these equations. 

The openings angle of the beam can be derived by ap- 
plying the cosine rule for sides (Eq.|A2j on the triangle plot- 
ted in the bottom panel of Fig. lAll 

cos p — cos a cos ^ -I- sin a sin cos ( A^) (A3) 

l|Gil et al.lll98i ). Here A0 is the pulse longitude measured 
from the fiducial plane. In the case of a profile which is 
symmetric around the fiducial plane A<j) = W/2. 

The azimuthal angle ip can be derived by using Eq. lAll 
on the same triangle plotted in the bottom panel of Fig. lAll 
and one directly obtains 

sin sin A0 



sm ifi - 



sm p 



(A4) 



l| Gupta et al.l |2004| ). It can be shown that this equation 
is equivalen t to th e following expression which is used by 
IWang et all l|2006h : 



tan 



sin (p — a) sin (p — p) 



sin p sin (p — Q 



(A5) 





Figure Al. Top plot: A spherical triangle. The angular lengths of 
the sides of the triangle are the angles a, b and c, while A, B and 
C are the angles between the sides of the triangle measured on 
the surface of the sphere. Bottom plot: A triangle on the celestial 
sphere used to derive the half openings angle of the beam p and 
the azimuthal angle of a fieldline ip. The line of sight makes an 
angle f with the rotation axis and and the magnetic axis makes 
an angle a with the rotation axis. 



where p = ^ {a + ( + p). 



